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Abstract 

Reaction and elastic differential cross sections are calculated for light nuclei in the framework 
of the Glauber theory. The optical phase-shift function is evaluated by Monte Carlo integration. 
This enables us to use the most accurate wave functions and calculate the phase-shift functions 
without approximation. Examples of proton nucleus (e.g. p- 6 He, p- 6 Li) and nucleus-nucleus (e.g. 
6 He— 12 C) scatterings illustrate the effectiveness of the method. This approach gives us a possibility 
of a more stringent analysis of the high-energy reactions of halo nuclei. 
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I. INTRODUCTION 



The study of hundreds of new unstable nuclei has become possible in the new radioactive 
beam facilities. The measurement of the cross sections at high energies is an important 
experimental tool for studying these elements PJ. The observed interaction cross sections 
can be related to the wave functions of these nuclei through the Glauber theory 0; thus 
one can obtain information on the structure of these exotic isotopes. This direct relation 
between the structure and reaction cross section is, however, hampered by several problems. 
The first and most serious one is that the calculation of the cross section by using accurate 
nuclear wave functions is difficult. Also the solution of the nuclear many-body problem is 
notoriously complicated. 

In Glauber theory the nucleus-nucleus elastic scattering amplitude is obtained by inte- 
grating the optical phase-shift function over the impact parameter. The optical phase-shift 
function is the matrix element of the multiple-scattering operator between the product of 
wave functions of the target and projectile. The difficulty of the evaluation of this matrix 
element stems from the fact that the multiple-scattering operator is an A-body (A is the 
total number of nucleons involved) operator, a product of pairwise nucleon-nucleon scatter- 
ing operators. To avoid the calculation of the matrix element of an A-body operator several 
simplifying approximations have been introduced over the years. One popular method is to 
replace the wave function of the target or both the projectile and the target by densities. 
This is then further simplified by assuming that the density of the nucleus is a product 
of one-body densities. In that case the matrix element of the multiple-scattering operator 
becomes a product of two-body matrix elements and is easy to calculate. Another popu- 
lar approach is to use an inert core and concentrate on the effect of the valence particles 
only. One can also expand the multiple-scattering operator as a sum of two-, three-, A- 
body operators and truncate the expansion at some level. The most popular method is the 
"optical-limit approximation" (OLA) where the complicated multiple-scattering operator is 
replaced by a simple two-body operator in a somewhat ad hoc way. 

The problem of these approximations is that as the "exact" scattering amplitude is not 
available, it is very difficult to judge how good they are. Some of the assumptions sound 
very serious; the usage of one-body densities most definitely washes out the effects of nuclear 
correlations. Some of these approximations may be especially bad for weakly bound halo 
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nuclei as has been discussed by many authors || f|, [5|, |6|, [7j . 

To avoid the complication of the analytical calculation of the cross section we will use 
Monte Carlo integration. We will show that the accuracy of this approach is very good 
and one can obtain very reliable results with very little effort compared to the previous 
approaches. Moreover, unlike the previous calculations, no approximation has to be in- 
troduced in the evaluation of the Glauber amplitude. We use many-body wave functions 
directly without loosing delicate correlation effects. In the framework of the Glauber the- 
ory, Monte Carlo integration has been used to calculate the nuclear transparency in the 
7n — > 7i~p process in 4 He and 16 in Ref. ||. 

This approach also allows us to use the most advanced quantum Monte Carlo (QMC) 
wave functions available for light nuclei. These wave functions are obtained by the solution 
of the nuclear Hamiltonian with realistic two- and three-body interactions. These ab initio 
QMC calculations of the energy spectra and other nuclear properties are in good agreement 
with experiments in the A < 8 region ||. The QMC predicts the radii of the light nuclei 
quite reliably so that a direct comparison of the nuclear sizes and reaction cross sections 
becomes possible. 

The road map of this paper is as follows. First we introduce the most important quantities 
of the Glauber theory in Sec. II. Next we show some of the details of the calculation in Sec. 
Ill and of the wave functions in Sec. IV. The presentation of the results can be found in 
Sec. V. A brief summary and outlook closes the paper in Sec. VI. 

II. GLAUBER THEORY 

Glauber's multiple-scattering theory of high-energy collisions, which is based on the 
eikonal and adiabatic approximations, is widely accepted as a standard tool for calculating 
various cross sections. The probability that two colliding nucleons in the projectile and 
the target lead to the excitation of the nuclei is almost unity in high-energy reactions (at 
several hundred MeV/A). Since the nuclear force is short-ranged, the probability of nuclear 
excitation in the collision reflects the geometrical size of the nuclei. The interaction or 
reaction cross section can thus be related to the size and structure of the nuclei. 

The basic quantity in the Glauber theory is the optical phase-shift function Xci(b) defined 
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by 



e*-<» = n nt 1 - + s r - i $ o (i) 



where 6 is the impact parameter, and $q ($q ) is the intrinsic projectile(target) wave function 
with its center-of-mass wave function removed. The profile function T is a two-dimensional 
Fourier transform of the nucleon-nucleon scattering amplitude 

/(M = £ fdbe-^ b T(b). (2) 
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In Eq. flU) one integrates over all (independent) intrinsic coordinates rf and rj. The two- 
dimensional vectors, sf and sj, are the projections of rf and rj onto the xy-plane which 
is perpendicular to the incident direction of the projectile. 

The nucleus-nucleus elastic scattering amplitude is easily obtained once the optical phase- 
shift function is available: 

/ cl (M) = ^ [dbe-^ b {l-e^% (3) 
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where K is the wave number of the relative motion between the two nuclei. The effects of 
the Coulomb interaction can be easily incorporated as well, but, as they are important only 
at extreme forward angles, we omitted the Coulombic contribution in the present analysis. 
The total reaction cross section is obtained by subtracting the elastic cross section from the 
total cross section: 



a r 



J db (l- |e %1 ^| 2 ). (4) 



The Glauber theory is a non-perturbative theory; its strength is that it directly employs 
information on the bare nucleon-nucleon interaction and thereby makes it possible to relate 
the reaction dynamics to the underlying interaction of the constituent particles. At lower 
energies and for heavier nuclei than those considered in this paper, the bare nucleon-nucleon 
cross sections must be corrected for medium effects related to Pauli blocking and effec- 
tive masses [f[0[| . Usually T is not calculated from the bare nucleon-nucleon interaction but 



parametrized in a convenient form (like Gaussians) so as to fit the empirical nucleon-nucleon 
scattering amplitude through Eq. (0): 



r( 6 ) = i ^-(l- ittT )expl--| (5) 



4 



where the parameters depend on the isospin of the nucleons (r=pp, nn, pn). Other operators, 
like spin-orbit and tensor, may also be necessary (especially if the spins of the target and 
projectile are nonzero), but those are most often neglected. Their inclusion would not cause 
any problem in the present approach. 

As seen in Eq. ([!]), the optical phase-shift function is defined by a many-body multiple- 
scattering operator and obviously its calculation is very involved. One often uses some 
simplification like the OLA, which is just the first-order approximation in the cumulant 
expansion [§]: 



e *xoLA(b) — ex p J _ / / dr p dr T p p (r p )p T (r T )T(s F - s T + b) } , (6) 



where e.g., p p (V p ) is the single-particle density of the projectile nucleus. 

Several authors have shown that a treatment beyond the OLA is necessary for a quantita- 
tive analysis of the reaction cross sections || f|, ^| as well as the elastic scattering cross sec- 
tions |], This is particularly true for loosely coupled systems such as halo nuclei because 
breakup effects are not properly accounted for in the OLA. There have been considerable 
efforts which attempt to calculate higher-order or all-order terms of the multiple-scattering 
operator. It is known that the matrix elements of higher-order terms should be evaluated 
by two- and more-particle densities 0. Using the density constructed from uncorrelated 
wave functions ]TT], [12], [13], |T3| is thus not sufficient to understand the importance of the 
correlated motion in nuclei even if all-order terms are evaluated. 

Very recently a method of calculating the complete Glauber amplitude |15| has been 



proposed and applied to p+ 6 He elastic scattering using a microscopic a+n+n three-cluster 
wave function for 6 He |[L6| . It uses an expansion of the multiple-scattering operator and 
evaluates each term analytically provided that the profile function T is expressed by Gaus- 
sians. The method has been extended to a nucleus-nucleus well with considerable 

success 



17]. However, this method is limited by two things: one is that the number of terms 
in the expansion becomes prohibitively large for heavier nuclei; another is that if T is not a 
Gaussian then the analytic integration for correlated wave functions is hopeless. 
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III. MONTE CARLO INTEGRATION 



Our purpose is to present a powerful method of calculating the matrix element of the 
multiple-scattering operator completely by taking another route, a Monte Carlo integration 
(MCI) of the multiple-scattering operator. No method has so far been available to calculate 
the Glauber amplitude completely for general correlated wave functions. The Monte Carlo 
integration seems to be the most natural way to calculate the phase-shift function. Its 
advantages are quite clear: (1) there is no restriction on the form of the target or projectile 
wave functions (general few- or many-body wave function can be used); (2) the full multiple- 
scattering operator can be used (there is no need for expansion or truncation); (3) it is very 
simple compared to earlier calculations. 

In our present approach the multiple-scattering operator is assumed to be independent of 
spins, depending only on the spatial coordinates and the isospins. The isospin dependence 
arises when different profile functions are used between the pp (nn) and pn pairs. The wave 
functions $q and $q are represented by a multicomponent (approximately 2 x (§)) vector 
in the spin-isospin space. The matrix elements are calculated by taking the scalar product 
of the bra and the ket in the spin-isospin space and using MCI in the coordinate space. 

In MCI, the integration points are generated by the Metropolis algorithm using |$q$q | 2 
as a guiding function. One may possibly get better accuracy and convergence by using 
(a+br 2 ) |$q $q | 2 (r is the root mean square radius of the projectile) or some similar expression 
which increases the weight of the asymptotic part, but as will be shown later the simple 
I^q^q | 2 form is sufficient for the present purposes. 

In the MCI we first generate a set of N integration points by a Metropolis random walk 
and store them. Then the optical phase-shift function is calculated over these configurations 
for each discretized value of the impact parameter b. In this way we not only save compu- 
tational time but we have a "correlated sampling" for different impact parameters avoiding 
the independent statistical errors of multiple Metropolis walks. The reaction cross section 
or the elastic differential cross section is calculated by numerical integration over the impact 
parameter. 

The impact parameter 6 is a two-dimensional vector which can be parametrized as 
(6sin0, 6cos0). For spherical nuclei the phase-shift function has no dependence on 0, so 
the integration in Eqs. ([3]) and (|4]) over b can be reduced to that over the radial variable b. 
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In the nonspherical case, one has a two-dimensional integration over the impact parameter, 
which increases the number of discretization points. The practical example studied in this 
paper shows that the dependence on is almost negligible. 

To test the MCI evaluation of the optical phase-shift function, a simple but nontrivial 
example, a+a scattering at 5.07 GeV/c (T Q =2.57 GeV), has been considered. Taking a single 
harmonic-oscillator shell-model wave function for the a-particle, the phase-shift function and 
the reaction cross section can be analytically calculated The MCI results are tested 

against this analytical example. Figure 1 shows the elastic differential cross section as a 
function of four-momentum transfer squared —t = h 2 q 2 using 100,000 MCI points. For small 
momentum transfer (this is the region where the present day experiments can be performed) 
a very small number of MCI points is sufficient to get reliable results. For larger momenta 
the oscillatory behavior coming from e~ iq b becomes more rapid, requiring more accurate 
integration for relatively large b, but even for the 10 5 points used the computational time is 
still almost negligibly small. Table I shows the reaction cross section and the mean-square 
(rms) radius of the a-particle obtained from this calculation. Naturally, the analytically 
calculated values are perfectly reproduced provided one uses a sufficiently large number of 
MCI configurations. These tests show that one can confidently use the MCI in calculation 
of the optical phase-shift function and related quantities. 



IV. WAVE FUNCTIONS 

The wave functions, except for 12 C, are obtained by the QMC method [^, [L8| using the 
Argonne vi& two-body ]19] and the Illinois three-body interaction IL2 p(J . The QMC method 



with these interactions provides a reliable description of the energy levels and different 
properties of light nuclei. Previous calculations are based on cluster, few-body, shell-model 
or mean-field wave functions to calculate the optical phase-shift functions. These calculations 
have used schematic effective interactions with adjustable parameters and other simplifying 
approximations. The use of QMC provides us with a realistic, ab initio, microscopic wave 
function. The 12 C nucleus is not yet accessible in QMC and for the 12 C a three-a microscopic 
cluster-model wave function is used. In this model the internal wave functions of the a- 
particles are single shell-model Slater determinants and the relative motion between the 
clusters is represented by linear combinations of Gaussians. The combination coefficients 
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are determined variationally by solving the 12-nucleon Schrodinger equation with an effective 
(Minnesota [pT| ) two-nucleon interaction. 

The QMC methods include variational Monte Carlo (VMC) and Green's function Monte 
Carlo (GFMC) methods. The VMC is an approximate method that is used as a starting 
point for the more accurate GFMC calculations. The VMC method starts with the con- 
struction of a variational trial function of specified angular momentum, parity and isospin, 
^T{J n ]T), using products of two- and three-body correlation operators acting on a fully 
antisymmetrized set of one-body basis states. Metropolis Monte Carlo integration is used to 
evaluate {9t\H\^t), giving an upper bound to the energy of the state. The GFMC method 
is a stochastic method that systematically improves on \I/<r by projecting out excited state 
contamination using the Euclidean propagation ^(t) = exp[— (H — E )t]^t- in the limit 
t — ► oo, this leads to the exact ground state energy. Details of the structure calculations 
can be found in ||. 

The calculation of expectation values directly using 

(Mt)\0\Mt)) 

( ( )} = <*M|*(t)> ' (?) 
is complicated and computationally too demanding. In our calculations we have used the 
approximate expression: 

(O(t)) « (0(r)> M ixed + [(0(r)) M i X cd - (0) T ] , (8) 
where the "mixed" expectation value between ty? and \P(t) is: 

(0)Mixed - (9 t \9(t)) ' (9) 

and (O)t is just the variational expectation value. This approximation is very good, provided 
that the difference between the VMC trial function ty? and the exact wave function is 
reasonably small. 



V. RESULTS 



In this section we present our results for total reaction and elastic differential cross sections 
of different projectile-target systems. The proton, neutron and matter rms radii of the nuclei 
predicted by the structure calculations are shown in Table II. The VMC gives a very good 
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wave function for the alpha particle, but underbinds the A=6 nuclei. The GFMC energies 
of 6 He and 6 Li are in very good agreement with experiment. The GFMC improvement of 
the wave function is especially important for the loosely bound 6 He. 

As a first application we have calculated the differential cross section for p+a elastic 
scattering at T p =0.7 GeV (see Fig. 2). (The parameters of the profile function (||) are taken 
from [T5|] ). The agreement with the observed cross section ]22|, |2Bj is perfect, but unfortu- 
nately the experimental data are only available up to 0.05 (GeV/c) 2 . The differential cross 
section obtained by a simple harmonic-oscillator a-particle wave function is also included for 
comparison. The realistic and simple shell-model prediction is nearly identical in the exper- 
imentally accessible region. For larger momentum transfer the two wave functions predict 
significantly different cross sections despite the fact that they give the same nuclear radius. 

Next we consider the a+a elastic scattering at 5.07 GeV/c. In this case the experimen- 
tal data |24| are available for a wider range of momentum transfer (see Fig. 3). While 
the simple shell-model wave function fails to explain the observed data points, the realistic 
wave function leads to good agreement with the experiments. The disagreement between 
the simple shell-model and realistic wave function prediction is almost an order of mag- 
nitude. To explain the data with a simple shell-model wave function one has to increase 
the harmonic-oscillator size parameter drastically and that would lead to an unrealistic a- 
particle radius (about 20% too large). This example clearly shows the importance of the 
realistic wave functions and the sensitivity of the experimental results to the details of the 
nuclear structure. 

Table III compares the reaction cross sections calculated for various systems with mea- 
sured interaction cross sections, a int . The interaction cross section is defined as the sum 
of the cross sections for all channels in which the nucleonic composition of the projectile 
changes. In high-energy collisions, the projectile can only lose nucleons, that is, the proba- 
bility of pickup processes is negligible, and the difference between cr rea c and <Ti nt is expected 
to be small. Possible differences between them come from two processes: one is the inelastic 
excitation of the projectile, which may occur if the projectile has a particle-bound excited 
state. Another is a process in which the projectile remains in its ground state while the tar- 
get gets excited. These processes contribute to cr reac but not to a int . When the projectile is 
a halo nucleus that has no particle-bound state, like 6 He, the second process can be ignored 
because the halo nucleus is easily broken by a small shock so it is unlikely that it remains in 
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its ground state while exciting the target. The example of 6 He+ 12 C collision at 0.8 GeV/A 
bears out this argument; the calculated reaction cross section is indeed close to the measured 
interaction cross section. The 6 Li+ 12 C reaction cross section is slightly larger than the inter- 
action cross section. It remains an open question whether this difference can be accounted 
for by the above processes. A similar comment may be applied to the 12 C+ 12 C case. For 
the sake of reference the reaction cross sections measured at 0.87 GeV/A are included in the 
table. 

Figures 4. a and 4.b present the elastic differential cross sections for p+ 6 He and p+ 6 Li at 
T p =0.7 GeV. The experimental data ^5j are available only up to —t = 0.05 (GeV/c) 2 . 
In that region (Fig. 5.b) the best theoretical (GFMC) cross sections slightly underestimate 
the experimental data, especially for 6 He. An ct+n+n cluster-model result [15| for p+ 6 He 
at T p =0.717 GeV, which is obtained by using the wave function [16| solved in a restricted 
model space with the Minnesota effective interaction, is also included for comparison. The 
cluster-model result agrees perfectly with the experimental data. (See Fig. 2 of Ref. 
The matter rms radius 2.51 fm in the cluster model, 2.56 fm in the VMC and 2.61 in the 
GFMC calculation. The slight underestimation of the differential cross section by the GFMC 
might be due to the fact that the size of 6 He given by GFMC is a little too large. For higher 
momentum transfer, there is a substantial difference between 6 He and 6 Li as well as the 
cluster and QMC results. 

The full and OLA calculations are compared for the reaction and elastic differential cross 
sections. Table III compares the reaction cross section for p+ 6 He and p+ 6 Li at T p =0.7 GeV. 
The OLA cross section is slightly smaller than that of the full calculation, which, differently 



from a nucleus- nucleus case, holds true for a proton-nucleus system ]|, pL5| . Figures 5. a and 
5.b compare the elastic differential cross sections of p+ 6 He and p+ 6 Li at T p =0.7 GeV. The 
difference in the cross sections at small four-momentum transfer between the full and OLA 
calculations is magnified in Fig. 5.b. Some difference can be seen in the cross section versus 
—t slope for the case of 6 He or in the magnitude of the cross section for the case of 6 Li at 
small \t\, but the more conspicuous differences appear near and beyond the first dip of the 
cross section. 

Differential cross sections versus four-momentum transfer squared are plotted in Fig. 6 for 
elastic scattering of several nuclei on a 12 C target at 0.8 GeV/ A incident energy. The angular 
distributions show Fraunhofer-type diffraction patterns and the first dips move to smaller 
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angles with increasing mass number. The basic feature of these trends can be understood 
in a strong absorption model (SAM), which is quite reasonable for high-energy collisions of 
nuclei. In an extreme version of the SAM, the phase-shift function is assumed to satisfy the 
relation 

• ™ f b < R 

e *xei(6) = 9(6 - E) = < , (10) 

I 1 b > R 

where R is a cut-off radius, and is on the order of the sum of the radii of the two colliding 
nuclei. (In reality e lXcl( ^ has a smooth cut-off distribution.) In the SAM the incoming flux 
corresponding to the collision with b < R is completely absorbed, while the collision with 
b > R receives no effect. (Note that the Coulomb interaction is ignored.) Substitution 
of Eq. (Top into Eq. (|]) leads to the scattering amplitude of f e \(9,4>) oc Ji(qR)/qR. The 
differential cross section \f e \{0, 4>)\ 2 thus vanishes at the zeros of the Bessel function J\(qR), 
whose zeros occur at qR ~ 3.83,7.02,.... With increasing A, R in general increases, thus 
the first zero appears at smaller values of —t. The value of R extracted from the dip in the 
figure corresponds to the empirical radius rather well. The 6 He and 12 C case does not follow 
this rule as the size of 6 He is larger than that of 12 C. This is mainly due to the fact that the 
sharp cut-off assumption is not good for 6 He because of its halo structure. The momentum 
transfer squared corresponding to the second dip is predicted as (7.02/3. 83) 2 ~ 3.4 times 
that of the first dip in the sharp cut-off model. 



VI. SUMMARY AND OUTLOOK 



The Monte Carlo integration was used to facilitate the evaluation of the complete Glauber 
amplitude involving the multiple-scattering operator. There was no need to introduce any 
ad hoc approximation or assumption in this approach. The great advantage of this method 
is that it enables us to use accurate, sophisticated wave functions of projectile and target 
nuclei. A number of calculations confirmed that it is possible to directly relate the wave 
functions to reaction cross sections measured at high energy. 

The calculations presented here focused on light nuclei. One can carry out similar inves- 
tigations for heavier elements provided that suitable wave functions are available. One very 
often replaces the wave function of the target nucleus with a density distribution. This ap- 
proximation renders calculations for heavier targets ( 27 A1 or 208 Pb) accessible. In that case 
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one has to construct a reliable nucleon-target profile function to replace the nucleon-nucleon 
profile function, but the rest of the calculation is the same as presented in this work. One 
further possible approximation is to use an inert core with correlated valence nucleons to 
construct the wave function of the projectile. Such a step might be necessary to investigate 
the cross sections of oxygen or sodium isotopes, for example. 

The experimental results mostly cover the low-momentum transfer region where the elas- 
tic differential cross section is not too sensitive to the details of the wave functions and simple 
models do quite well. We have shown, however, that there is a strong discrepancy between 
the predictions of simple model and realistic wave functions. The high-energy reactions 
therefore may give information on the details of the nuclear structure. For example, the 
a + a scattering at 5.07 GeV/c, where the experimental data is available for larger momen- 
tum transfers, can only be described by using realistic wave functions for the alpha particle. 
One hopes that the new radioactive beam facilities will provide us experimental data at 
higher momentum transfer and further tests of the wave functions will become possible. 
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TABLE I: Comparison of a + a reaction cross section at 5.07 GeV/c and the rms radius of the 
a-particle calculated analytically and by Monte Carlo integration with 100,000 points. The wave 
function of the a-particle is taken as the simple (0s) 4 Slater determinant. The statistical error of 
the Monte Carlo integration is given in parenthesis. 

Method cr rcac (mb) (r 2 ) (fm 2 ) 

analytic 242.91 2.250 

MCI 242.8(7) 2.251(3) 

TABLE II: Point proton, neutron and matter rms radii and reaction cross sections for the collision 
with a 0.7 GeV proton. The Monte Carlo statistical errors for <7 reac are 1 mb in all cases. 



System 


wave func" 


tion r p (fm) 


r n (fm) 


r m (i 


in) o- rca , 


a 


VMC 


1.46 


1.46 


1.46 


100 


6 He 


VMC 


1.96 


2.81 


2.56 


163 


6 He 


GFMC 


1.92 


2.87 


2.61 


172 


6 Li 


VMC 


2.47 


2.47 


2.47 


165 


6 Li 


GFMC 


2.47 


2.47 


2.47 


166 


12 C 


cluster 


2.36 


2.36 


2.36 


254 
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TABLE III: Calculated reaction cross sections. Experimental data are the interaction cross sections 
taken from [26|. Cross sections with t are the reaction cross sections at 0.87 GeV/yj f2^| . The 
statistical error of the Monte Carlo integration is given in parenthesis. The reaction cross section 
marked by * are obtained in the OLA case. 



bystem 


Energy (GeV/A) 


<7 rea c(mb) (Tint (mb) (Exp.) 


p+a 


0.7 


100(1) 




p+ 6 He 


0.7 


172(1) 




p+ 6 Li 


0.7 


166(1) 




p+ 6 He 


0.7 


166* 




, fir - 

p+ b Li 


0.7 


164* 




^He+^C 


0.7 


721(2) 




b Li+ 12 C 


0.7 


lUo{z, ) 




a+a 


0.6425 


235(1) 




P+ 12 C 


0.8 


261(1) 


262±13.5 f 


a+ 12 C 


0.8 


506(1) 


503±5, 527±26t 


p+ 6 He 


0.8 


178(1) 




p+ 6 Li 


0.8 


171(1) 




6 He+ 12 C 


0.8 


733(2) 


722±6 


6 Li+ 12 C 


0.8 


712(2) 


688±10 


12 C+ 12 C 


0.8 


865(1) 


856±9, 939±49 t 
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FIG. 1: Center-of-mass differential cross section versus four- momentum transfer squared for a + a 
elastic scattering at 5.07 GeV/c. The solid curve is the result of the analytical calculation, while 
the result of Monte Carlo integration is shown for 100,000 points. 
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FIG. 2: Comparison of the differential cross section for p+a elastic scattering at T p =0.7 GeV 
obtained by a simple shell-model (dashed line) and a realistic (solid line) wave function of the 
a-particle. The two wave functions give the same rms radius; Monte Carlo error estimates are 



shown for the latter case. The data are available only for -t < 0.05 (GeV/c) 2 ||, ||, but not 
shown because they are hardly distinguishable from the theoretical curves. 
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FIG. 3: Comparison of the experimental and theoretical differential cross section for a + a elastic 
scattering at 5.07 GeV/c. The solid line shows the results with a realistic a-particle wave function, 
while the dashed line is obtained by a simple 0s shell-model wave function. The two wave functions 



give the same rms radius. The data are taken from [24| 
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FIG. 4: (a) Differential cross section versus four-momentum transfer squared for p+ 6 He and p+ 6 Li 
elastic scatterings at T p =0.7 GeV. An a+n+n cluster-model result (T^j is shown for p+ 6 He at 
T p =0.717 GeV. The statistical error of the Monte Carlo integration is shown for the GFMC and 



VMC results, (b) Theoretical and experimental [25] differential cross sections for p+ 6 He and p+ 6 Li 
elastic scatterings at T p =0.7 GeV. 
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FIG. 5: Comparision of the full and OLA differential cross sections. The parameters are listed in 
the caption of Fig. 4. The GFMC results are used to calculate the cross sections. 




FIG. 6: Differential cross sections versus four-momentum transfer squared for elastic scatterings 
of p, a, 6 He and 12 C on 12 C at 0.8 GeV/A 
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